Journal of Nonlinear Mathematical Physics 1998, V.5, N 1, fllH53[ Article 

Progressive Internal Gravity Waves With 
Bounded Upper Surface Climbing a Triangular 
Obstacle 

Mina B. ABD-EL-MALEK ] and Malak N. MAKAR % 

f Department of Engineering Mathematics, Faculty of Engineering, 
Alexandria University, Alexandria 21544, Egypt 
E-mail: mina@alex.eun. eg 

\ Department of Mathematics, Faculty of Science, 
Assiut University, Assiut 71516, Egypt 

Received January 10, 1997; Accepted July 10, 1997 

Abstract 

In this paper we discuss a theoretical model for the interfacial profiles of progressive 
non- linear waves which result from introducing a triangular obstacle, of finite height, 
attached to the bottom below the flow of a stratified, ideal, two layer fluid, bounded 
from above by a rigid boundary. The derived equations are solved by using a nonlinear 
perturbation method. The dependence of the interfacial profile on the triangular 
obstacle size, as well as its dependence on some flow parameters, such as the ratios of 
depths and densities of the two fluids, have been studied. 

1 Introduction 

The determination of flow patterns over obstacles is a problem of a great interest, that 
attracted many scientists over the past decades. Lamb [6] was the first to give the essential 
features of the flow of an ideal fluid in an open channel in the presence of an obstruction 
in the channel. In 1955, Long [7] and then later on Mclntyre [8] in 1972, considered the 
case of a steady and uniform stratification over obstacles of finite height, while Mei and 
LeMehaute [9] in 1966 studied the case of long waves in shallow water over an uneven 
bottom. The effect of the irregularities of the bottom, on gravity waves, has been studied 
by Kakutani [4] in 1971 via a reductive perturbation method. Recently, Kevorkian and 
Yu [5], in 1989, studied the behaviour of shallow water waves excited by a small amplitude 
bottom disturbance in the presence of a uniform incoming flow. In this paper we study 
a theoretical model for the interfacial profiles of progressive non-linear waves result from 
introducing a triangular obstacle, attached to the bottom below the flow of a stratified, 
ideal, two layer fluid, bounded from above by a rigid boundary. Our primary motivation 



Copyright ©1998 by M.B. Abd-el-Malek 
and M.N. Makar 



42 



M.B. Abd-el-Malek and M.N. Makar 



for the present investigation is calculate the shape of the interfacial waves, and to discuss 
the influence of both geometrical and flow parameters on the profiles. In section 2, we 
extended the mathematical technique applied by Helal h Molines [3], 1981, in determining 
the nonlinear free-surface and interfacial waves in a tank with flat horizontal bottom. 
Nonlinear perturbation technique is used, leading, in sections 3 and 4, to expression for 
interfacial wave, was derived in the form of expansions in powers of where £ is a small 
parameter that provides a measure of weakness of dispersion. Boutros, et al [1], in 1991, 
applied the same technique to study the internal waves over a ramp. 

Finally, in section 5 we have discussed the effect of the density ratio, R, the thickness 
ratio, H, and the triangle height, L. 

2 Formulation of the problem 

Two-dimensional irrotational motion of a stably stratified two-layer fluid with a rigid 
upper boundary and a bottom surface in the form of a triangle with two inclination angles 
a and (3. A cartesian coordinate system is defined with the origin at the bottom surface. 

We assume that the motion is two-dimensional, and the fluid is inviscid, incompressible, 
and that the flow field due to the wave motion remains irrotational and consequently we 
can introduce velocity potentials of upper and lower layers are denoted by i = 1,2, 

respectively. Moreover, let H^, p^*); i = 1,2 denote the thickness and densities of the upper 
and lower fluids, respectively, r* is the time, Y* = W*{X*) is the bed of the channel and 
Y* = f*(X*,T*) is the interfacial disturbance from uniform condition. The component of 
gravity, vertically downwards, is g. The equations of motion are thus the Euler equations 
together with the continuity equation. All variables are non-dimensionalized by using the 
characteristic length i^l and time and accordingly 

U = U*/[gH^f'^ and = / {Hl[gH^Yl^') . 

Moreover, assuming that the fluids are in the undisturbed uniform state up/down stream 
at infinity, we impose the following boundary conditions with respect to X* 

^*^} = U*, (i = l,2) as X*^±oo. 

An essential step which makes our problem easier in handling is to define an appropriate 

stretching of the horizontal coordinate while leaving the vertical coordinate unchanged due 
to the fact that the horizontal dimensions are much greater than the vertical dimensions, 
thus we define 

X = eX, y = Y, t = £T, (1) 
where £ is a small parameter. Thus the basic equations for this system can be written as 
£2$« + $g=0, f<y<l + H, -^<x<^, (2) 

£2^2) + $(2) = 0, M^<y</, -oo<x<oo, (3) 
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with conditions 

(i) Boundary conditions: 



a = 1,2) 



.(1) 



1 



(2) , 1 

2 



2 



(2)^2 



at y = f 



^i^)=£^^pw^ at y = W{x) 
$« = 0, at y = l + H, 

= 1, {i = 1, 2) as x ^ ±oo 



(4) 

(5) 
(6) 
(7) 



(ii) Initial condition: 

at t = 0: the initial profile of the interfacial wave, denoted by f{x, 0), is shown in Fig.l. 

fix,0) 



1.005 



0.995 
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13, Xe = 30, L = 0.25, Q = 0.01923, 



Fig.l. The inital waveform over a triangular obstacle with x 
P = -0.01471 



where the density ratio R = p^^^ /p^^^ (less than unity) and the thickness ratio H are two 
characteristic parameters of the system, and W{x) has the form 

W{x) = ax + b, (8) 

where 

'(0,0) x<0 
(a, 0) < X < Xm 

( l^^e) ^ ^ ^ 

(0,0) X>Xe 

Since we consider weakly nonlinear waves, we expand the dependent variables as power 
series in the same parameter e around the undisturbed uniform state, following Helal and 
Molines [3], we get 



(a, 6) 



<^i.i) = Y^e''--^G^l_^{x,y,t), i = l,2 

n=0 
oo 

/ = E^'"/2n(x,y,t), with /o = l. 



n=0 
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The scale parameter e, which is assumed to be smah, provides a measure of weakness 
of dispersion. 

The boundary conditions on the interface, equations (^), are expanded as a Taylor 
expansion of the type 



n=0 



QyU 



(9) 



yo 



When dlj), (P), using the expansion (^), are inserted into equations (§)-(0) and powers of 
e are sorted out, we get an ordered set of equations to be solved. 



3 Orders of approximations 

3.1 The first-order approximation 

Equations of the first-order approximation, finally gives, for i = 1,2 

where B^'^\x,t) are unknown functions to be determined. 

3.2 The second-order approximation 

From the equations obtained from the second-order approximation, we conclude that 
5^=0, (« = 1,2) as x^±oo 

and 



f2ix,t) 



1 



l-R 



RE 



(1) 



B 



3.3 The third- and fourth-order appoximations 

Equations of the third- and fourth-order approximation, finally gives, for i = 1, 2 
Gf = -ly^B^l+yC^\x,t) + D^^{x,t), 



(10) 



where C^*^(x,t) and D^^\x,t) are arbitrary functions satisfy the following boundary con- 
ditions: 



C^*)=0 (i = l,2) as x^iboo, 
C7(2)(a;,t) = (WB^^ at y = W{x) 
dW=0 a = 1,2) as x^±oo. 



(11) 
(12) 
(13) 



Substituting equation (10) in the equations that are obtained from the third- and 
fourth-order approximation, we obtain 



(// + l)i?«-C«=0, 



(14) 
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and for i = 1,2 



B'i - C« + ^ (rB^^ - 4f ) = 0. (15) 



From equations (|12|), (Ml), and (|15D we get 

(16) 

where Qi, ^2 are the differential operators defined by 

U,^-H{l-R)—^+R—^, 
ox Ol 

52 dW d 

From equations (|l6|)-(p!7|) we can get, after getting rid of B^^^ and substituting for W{x), 
the following differential equation for the unknown function B^'^\x,t) 

-H{1 - R){1 - b - ax)Bi%, + [H + R{l-b- ax)]41 

-RaB^^^t + 3Ha{l - R)Bf^^ = 
and for /^(xjt) we can get the following relation 



(18) 



h{x,t) = -^\R 



3.4 The fifth- and sixth-order approximations 

Equations of the fifth- and sixth-order approximation lead to, for i = 1,2 

4'^ = 1^41. - jC^^i^^t) - ylD^il{x,t) +yE^){x,t) + F^^{x,t), (19) 

where E^^\x,t) and F^'^\x,t) are arbitrary functions, satisfy the folowing conditions: 

Ei!^=0 ii = 1,2), as x^iboo (20) 
and at y = W{x) 

F«=0 (i = l,2) as x^iboo. (22) 

Introducing equations (p^)-(p^) in the boundary conditions, we have the foUowind rela- 
tions: 

- ^-^^CW -iH + 1)41) + E^^^ = (23) 
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and for i = 1,2 



1 



1 



1 



(1)^ 5» 



2) 



(24) 



-R [ -WxL + ciP + A?^ ) + 4?^ - RB'J'B^,, 



.(1) , n(i) 



0. 



Thus the problem is now reduced to solving equations (14) and (|l^) for iJ^*-* and C*-*^ and 
next equations (||), (^) and (|§ for and where i = 1,2. 



4 Case of progressive wave 

It must be remarked that our procedure is valid as long as a ^ e^, otherwise a two- 
parameter analysis has to be carried out. Moreover, we shall invoke the smallness of a 
and write perturbation expansions for B^^\ {i = 1,2), in the form 



(*) _L ^2 



(25) 



Substituting (|25| ) in ([l8| ) and equating coefficients of a'^^\ j = 0, 1, 2, . . . we get the fol- 
lowing system of differential equations 



□5?) = AB^, 



(i = 0,1,2,...), i?i'}=0, 
where □, A are two differential operators defined as 



(26) 



Q3 



Equation (p^), for j = 0, has the following general solution, for the case of pure progressive 
waves, 



B, 



(0 _ R(i) 





with 



C = x--ft, 



^_H{l-b){l-R) 



H + {l-b)R ■ 
From equations (^), (p^), and ( p5[ ) we get 

oo 



n=0 



si^) = 0. 



(27) 



Again substituting equations (p5[), ( [27| ) in equation (|T5|) we get, after equating coefficients 
of a", n = 0,1,2,.... 



(2) 

0,x 



(2) 



BqI + A-BH, 



1-5 
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where 

The ehmination of E^^^ in equations (|2^ ) and (|2^ ) gives, for "a", the following system 
of differential equations 



^'R n(i) 



.(2) 



(28) 



1 - 



7' 



(1) 



(1) h(1) 



where 



Pi 



-H{2H^ + 6H + 3) ^ 7^[A( 1 - 26) + /?(2F + 1)] 



6 



A 



P2 = -(1 - 36 + 26^) + , 

6^ ^ 2(l-i?) 



2(1 -i?) 
[(26- l)A-i?(2/7 + 1)], 



Qi 



7 



(A2 + 2A-3i2), Q2 



1-R 



[i?(2A + l)-3A2]. 



For the non-trivial solution of D^^^ and D^^^ , the following differential equation for i?g^^ 
should be satisfied: 



(29) 



where 



Ml 



1-6 



7 



1-R 



Pi 



l-R 



-Pi, Mo 



1-6 



7 



1-i? 



Qi 



1 



Q2. 



Define 



r - 5^^) 



(30) 



Thus equation (|29D , by virtue of equation (|30[), will be transformed to the Boussinesq 
equation 



MiFg^g + MaFT^ = 0. 



(31) 



Helal & Molines [3] mentioned that the general solution of equation (^Tj) was found by 
Byrd and Friedmann [2] to be, in terms of the Jacobi elliptic function sn(u, /c), as 



0(1) _ V 



1 



3A;2 



where Yi is the greatest of the roots of the polynomial resulting from integrating equation 
(31) twice and k is the modulus of the Jacobean elliptic function, and 

1 / 3AYi \ 
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For small values of k the above elliptic function could be calculated in terms of trigono- 
metric functions, see Milne-Thomson [11], thus we have 



B, 



(1) 

0,? 



k 



1 



3k' 



A;2 + l 



k 



1 k' 

k^ 



16 



+ 



cos zo^ — cos 4o4 



128 



k^ 



— cosQ5i-5i{ ( — + — 1 sin2(5C + — sin450 + 5^^^ <! — + — cos 25^ 



64 

k^ k^ 



128 



16 



(32) 



(2') (2\ 

Substituting in equation (^) , for ^ and Bq / , we get the following fourth-order linear 
partial differential equation 



□^(2) = (AnX sin 2n(5^ + An+e cos 2n6£,) 



n=l 



+ (5^ ^ {An+3X COS 2n6S, + An+io sin 2n(5^) 



(33) 



n=l 



r {Aex sin 25^ ylia cos 26^) + Aio , 

where the coefficients Ai, A2, . . . , A13 are given at the end of the paper, as Appendix 1. 

Solving equation ( |33|) for the unknown B^\ following Miller [10], and calculating B^^ 
we get 



B 



(2) 



(34) 



Bq'^ + ^1*^ + r2x''t + (ra -|- r4x2 + rsxt + rgt^) sin 2(5^ 
_)_ (^fj + rgx^ -|- rgxi -|- not^) sin 45^ -|- rn sin 65^ 
+ (ri2 + ri3X + ri4t + ri^x^ + rigx^t -|- ri7a;t2 + rigt^) cos 2(5^ 
-I- {rigx + r2ot) cos 45^ + {r2ix + r22t) cos 66^, 

where the coefficients ri, r2, . . . , r22 are also given at the end of the paper, as Appendix 2. 
Taking into consideration the value of Bq^I from equation (^2|), we can get Bq'I and 



thus, using (^4|) for B\ [ we can get B\^l 



B 



(1) 



1 - a 



-B, 



(2) 
0,t 



In order to account for the nonlinear effects the 0{e'^) equations have to be considered 
as well. Thus bearing in mind the linear system of equations (28), the principal and 
secondary determinants of this system, we come to the result that 

dP=0, (i = l,2). 

Hence f4,{x,t) may be rewritten in the simplified form 

1 



f4{x,t) 



2(1 -R) 



{X-R) + 2(1 + H)R - 2X{ax + b) + 



ax\{l - 2b) 



1-6 



B, 



(1) 

O.xxt 



+{Ri2H + 1) + A(l - 2b))B[%, + {2{R- X') 



2ax\^{ax + 2)' 



(1)^2 



(35) 



+a\R-a){B^^lf + 



2aXb 



1 



B^% + 2a{R-X' 



axX 



-°0,x-°l,x 
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Hence f{x,t) will take the form 

where /4(x,t) is given by (^) and B^^ and B^^^ are given by 

5 Presentation of results and discussion 

The number of terms which has been obtained seems to be a good measure for the purpose 
of illustrating the effect of the parameters the density ratio, R, the thickness ratio, H, 
and the obstacle height, L. The error, difference between the fourth and second order 
approximations, in the interfacial profile for the two approximations is of order 10~^. Thus 
we limit our calculations up to the second-order approximation, as well as we considered 
the following values for the description of the triangular obstacle: Xm = 13 and Xg = 30. 

We studied the effect of the density ratio, R, on the wave profiles at the interfacial 
surface. Three values of R have been considered, namely R = 0.7, 0.8, and 0.9 for fixed 
values of H, L, and t. It is clear that as R decreases, there is a kind of violent oscillations in 
the obstacle region. This phenomena vanishes gradually as "i?" increases. An important 
remark must be mentioned is that, for the interfacial wave in the downstream region, the 
period of oscillation is much longer for the case when the two fluids are of very nearly 
equal density than that of significant different densities. This is due to the fact that the 
presence of the upper fluid has the effect of decreasing the velocity of propagation of the 
wave which consequently causes the decrease of the potential energy of a given deformation 
of the interface as well as the increase of the inertia. This result comes in good agreement 
with Lamb [6], who gave a marvelous natural example for such a phenomena, occurring 
near the mouths of some of the Norwegian fiord, when there is a layer of fresh water over 
salt water. 

The interfacial wave profiles, f{x,t), has been studied for different values of the thick- 
ness ratio, H, namely H = 0.3, 0.5, and 0.6 while the other parameters R, L, and t are 
fixed. It is clear that as H increases, there is an increase in the amplitude of the wave 
along the obstacle interval, as well as an increase in the wave length. 

We study the effect of changing the triangle height, L. Three values of L have been 
considered, namely L = 0.1, 0.2, and 0.25 for fixed values of R, H, and t. For the 
interfacial wave, as L increases a kind of violent disturbance in the wave profile appears, 
starting by a sudden increase in the profile, ending by a steep decrease at the beginning 
of the downstream interval. The behaviour of that solution can be interpreted, following 
Kakutani [4] , as follows: a given smooth waveform will propagate along the characteristic 
curves, gradually steepen its shape due to nonlinear interactions, and then the dispersive 
term will begin to play its role to balance this steeping. 

Appendix 1 

Ai = Wi (^-4 + 6^ + ^k^^ A2 = Wi{2k^ - 8k^) A3 = Wi ("^fc') 
= WiiAk"^ - 2k^) A^ = 4Wik^ Ae = Wik^ 



32) and (pi) respectively. 
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2 / ^Ylk^5^ 



where 

Wi = [7^ - H{1 - R)] 

and 



fc2 + l 



Ar = W2(2-2e- As = W2 (2k'' - ^fc^) = ^W2k^ 

AiQ = \w2k^ All = f 2A;2 - ^A;^^ A12 = W2k^ A13 = -^WaA;^ 



where 



and 



2 



Ai^ = H{l-R){b-l) Ai5 = H + l-b A16 = -^{2Ai - A^ - Ae) 

11 1 
An = ^SAi2 A18 = ^{Au - ^13) ^19 = Y^(4^2 - A5) 

A20 = A21 = -^As A22 = 7^(^13 - 2A7 - All) 

47 670 470 

A23 = -^(^4 + Aq) A24 = -^<^^13 ^25 = -Y^*-^"^^ ^ "^^^-^ 

47 D70 A15 — 27^ 7 
^30 = —2^28 ^31 = ^15^28 ^32 = ^15^28 ^33 = 7^15^28 

7 

^34 = 7^—^14^28 ^35 = 37^28 ^36 = 7^—^15^28 ^37 = 7^—^28 
27 27 27 

^38 = (7^^i5 + 2^14)^28 ^39 = ^{^^u + 7^^15)^28 A40 = 2^15^28 

27 

A41 = 7^(2^1 -A4-Ae + 2S[Aii - A^]) A42 = SAq A43 = \^5Aq 

470 z 

Am = l{Ai3 - All) A45 = -^(4^2 -^5 + 4Mi2) A^e = —A^ 
2 ID70 4 

A47 = -7^(^4+^6+5^13) ^48 = ^(^4 + ^6 + 25^13) ^49 = -j75^13 
Z7 z z 
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^50 = ^^5 -451 = SAls^as = 3^18^40 As = 6^8(^0-439 + A32A40) 

-454 = 6^8-430-435 -455 = 6^38 (-430 Ao + -432-435) -456 = -3^15^8 
A57 = -SAssiAlo + 2^39^5) ^58 = 2^35^40-438 

-459 = -429 + 2^31^8 - 2^9-440 + 3^38^2 + -457 -460 = -430 - -415 
Aei = A31 + 2^8-432 - -4lo - 2^39-435 + 3-418^30 + ^58 
-462 = -432 + 2^8-430 - 2^35^40 + -456 

^63 = 2(^38^36 + -432^0 + -432-439) - -440 - 6A39A35 + -453 

-464 = -436 + 2(-43i^35 + ^32-44o) - 3^5 (-435-439 + ^lo) + ^55 

^65 = 2(yl30-440 + -432-435) + -454 - 3^15^0 -466 = -439 + 2^8-440 + -451 

-467 = -440 + 2^38-435 -468 = 2^38-439 + -452 

-469 = -425 + (2^1^126 + -462-450) -470 = -450 + 2^38-426 

-471 = ^-438 (-450 + 2^8-426) -472 = ^(-435-450 + 2^26-467) 

-473 = ^ (2-466-426 + -467-450) -474 = -422 + 2^ (-460-449 + -461-447 + 2^62-448) 

-475 = -448 + 2^8-447 -476 = -449 + ^-438 (-448 + 2^38^47) 

-477 = ^ (-435-448 + 2^7-447) -478 = ^ (-435-449 + 2^67-448 + -466-447) 

-479 = ^(-435-446 + -467-445) -480 = -441 + ^(^0-443 + 3^61-417 + -462-442) 
-481 = -442 + 3^8-417 -482 = -443 + -438-442 + 3-4|8-4i7 
-483 = -444 + -438-441 + ^(-461-442 + -462-443) 

-484 = ^(-438-443 + -418-442) + -4i8-4l7 

-485 = ^(3^63^17 + ^64-442 + -465-443) + ^(-435-444 + -467^0 

^86 = ^ (3^66-417 + -435-443 + -467-442) -487 = ^ (-466-442 + -467-443 + 3^68-417) 
-488 = 777 (-435-442 + 3^67-417) 
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Appendix 2 

ri = -^^10 r2 = ^A37[2-fS{Au - ^ss) - ^78 - ^ss] 

rs = g^^37[27<^(-447 - ^88) - -48i] r4 = ^^37[7'^(^75 - ^86) - ^82] 
r-5 = ^^37[27<^(^76 - ^87) - 3^84] 

re = ^^^37 [47(^(^69 - ^79) - ^73 - ^46 - ^38^45] = ^^7^37^26 

'^8 = ^^7^37^0 rg = Y^7^37^71 

riQ = 216(53 ^37(7(6^27 - ^67) - ^38^2i] = ^^37^80 

ri2 = ^37(27(^^77 + A75 - Ase) ri3 = ^^37(7^(^73 + ^83) + ^76 - ^87] 

ru = ^7-437-488 ri5 = ^7^37^181 ne = ^7^37^182 

ri7 = ^7^37-484 ri8 = ^^^37(^70 + 47(5(^72 + A45)] 

ri9 = ^^^37[^71 + 27(5(^73 + ^46 + ^38-445)] ^20 = ^^7^37^21 
r21 = ^^7^37^38^21 r-22 = ^^7-438-421^43 
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